Finite electron crystallites in strong magnetic fields: Precursors of a supersolid? 
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We show that a supersolid phase, exhibiting simultaneously solid and superfluid behavior, properly 
describes the finite electron crystallites that form in two-dimensional quantum dots under high 
magnetic fields. These crystallites rotate already in their ground state and exhibit a nonclassical 
rotational inertia. They are precursors to a supersolid crystal in the lowest Landau level. We use 
exact numerical diagonalization, calculations employing analytic many-body wave functions, and 
a newly derived analytic expression for the total energies that permits calculations for arbitrary 
number of electrons. 
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The existence of an exotic supersolid crystalline phase 
with combined solid and superfluid characteristics has 
been long conjectured This conjecture was made 

in relation to solid 4 He. The recent experimental discov- 
ery Q that solid 4 He exhibits a nonclassical (nonrigid) 
rotational inertia (NCRI has revived an intense in- 
terest [E IE IE E| m the existence and properties of the 
supersolid phase. 

The supersolid phase is not restricted to bosonic crys- 
tals 9, 101. In this paper, we demonstrate the existence 
of such a phase in two-dimensional (2D) semiconductor 
iV-electron quantum dots (QDs) under strong perpen- 
dicular magnetic fields (B). Central to the concept of a 
supersolid is the appearance of a NCRI when the crystal 
is made to rotate |E IE] • 

Under a high magnetic field, the electrons confined 
in a QD localize at the vertices of concentric polygo- 
nal ring s and form a rotating electron crystallite (REC) 
[ill |l2j . The ground-state many-electron wave function 
has a nonzero good total angular momentum quantum 
number. We show that the corresponding rotational in- 
ertia strongly deviates from the rigid classical value, a 
fact that endows the REC with supersolid characteris- 
tics. The REC at high B can be viewed as the precursor 
of a supersolid crystal that develops in the lowest Landau 
level (LLL) in the thermodynamic limit. These conclu- 
sions were enabled by the development of an analytic 
expression for the energy of the REC that permits cal- 
culations for an arbitrary number of electron s, g iven the 
classical polygonal-ring structure in the QD [13j . 

For strong B we can approximate |l4| the single- 
electron wave function by (parameter free) displaced 
Gaussian functions; namely, for an electron localized at 
Rj (Zj), we use [l4( the expression 
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with A = y/ti/m*Cl; = ya 



?/4, where 



eB/(m*c) is the cyclotron frequency and uoq specifies 
the external parabolic confinement defining the QD. 
The position variables are given by the complex num- 



ber z — x + iy and Zj = Xj + i Yj . The phase guaran- 
tees gauge invariance in the presence of a perpendicular 
magnetic field and is given in the symmetric gauge by 
(p(z,Zj-,B) = (xYj - yXj)/2l%, with l B = y/hc/eB be- 
ing the magnetic length. 

For an extended 2D system, the Zj's form a triangu- 
lar latti ce Il5l. For finite N, however, the Zj's coin- 
cide 0, Il4 Il6| with the equilibrium positions [forming 
r concentric regular polygons denoted as (rii,7i2, ...,n r )] 
of N = Y^,q=i n q classical point charges inside an exter- 
nal parabolic confinement |13| . The wave function of the 
static electron crystallite (SEC) is a single Slater deter- 
minant |5 ,SEC [z]) made out of the single-electron wave 
functions u(zi, Zj), i = 1, N. 

| tj, sec j-^j ^ re p resen ts a broken-symmetry state, since it 
is not an eigenstate of the total angular momentum op- 
erator with eigenvalue L. For a rotating electron crystal- 
lite, L must be a good quantum number. To decsribe the 
rotation of the electron crystallite, we project the Slater 
determinant onto a state with good L 0, Q, 0, 0] , 
thus restoring the circular symmetry and obtaining 
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x|* SEC ( 7l ,..., 7r ))exp J . (2) 

Here L = Y^ q =i L q ancl | V I /SEC [7]) is the original Slater 
determinant with all the single- electron wave functions 
of the qth ring rotated (collectively, i.e., coherently) by 
the same azimuthal angle j q . Note that Eq. @ can 
be written as a product of projection operators acting 
on the original Slater determinant [i.e., on |5 ,SEC (7i = 
0,...,7 r = 0))]. Setting A = IbV% restricts the single- 
electron wave function in Eq. Q to be entirely in the 
lowest Landau level 11]. The continuous-configuration- 
interaction form of the projected wave functions [i.e., 
the linear superposition of dctcrmimants in Eq. (0)] im- 
plies a highly entangled state. We require here that B 
is sufficiently strong so that all the electrons are spin- 
polarized and that the ground-state angular momentum 
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FIG. 1: Ground-state energies calculated with the REC wave 
function compared to EXD calculations. Energies for N = 4 
fully polarized electrons (referenced to 4hQ) as a function 
of the magnetic field B. Dashed line (online red): Broken- 
symmetry SEC. Solid line (online green): EXD (from Ref. 
[lg). Dashed-dotted line (online blue): Symmetry-restored 
REC. Numbers under the curve denote the value of magic 
angular momenta L m of the ground state. Corresponding 
fractional filling factors are specified by v — N(N — l)/(2L m ). 
Parameters used: confinement hu>o — 3.60 meV, dielectric 
constant k = 13.1, effective mass m* — 0.067m e . 

L > La = N(N — l)/2 [the minimum value L$ specifies 
the socalled maximum density droplet]. 

The energy of the REC state [Eq. @] is given by [3, 

E3 

ET C = ^ 2 \([ 7 ]) e lW - [L1 d[ 7 ]/£^([7])e 4H - [il d[7], 

(3) 

with h([i\) = (* SEC ([0])|i/|* SEC ([ 7 ])), n([ 7 ]) = 
(* SEC ([0])|* SEC ([ 7 ])), and [ 7 ] • [L] = EUT-rAr The 
SEC energies are simply given by £sec — h([0]) / n([0]) . 
The many-body Hamiltonian is 

N i „ AT N N 9 

i— 1 i—1 i—1 j>i J 

(4) 

which describes N electrons (interacting via a Coulomb 
repulsion) confined by a parabolic potential of frequency 
luq and subjected to a perpendicular magnetic field B, 
whose vector potential is given in the symmetric gauge 
by A(r) = |(— By, Bx, 0). m* is the effective electron 
mass, k is the dielectric constant of the semiconductor 
material, and = |r< — rj|. For sufficiently high mag- 
netic fields, the electrons are fully spin-polarized and the 
Zeeman term (not shown here) does not need to be con- 
sidered. 

Prior to demonstrating that the REC wave function 
[Eq. (0)] exhibits supersolid characteristics, we address 
the numerical accuracy of its energies [E^ EC in Eq. J2J] 
compared to exact diagonalization (EXD) solutions of the 
many-body problem defined by the Hamiltonian in Eq. 
(0J. In Fig. 1, the ground-state energies as a function of 



B are compared to EXD calculations 18] for a QD with 
N = 4 electrons. The thick dotted line (online red) repre- 
sents the broken-symmetry SEC approximation, varying 
smoothly with B and lying above the EXD curve [solid 
(online green)]. The REC curve [dashed-dotted (online 
blue)] agrees well with the EXD one in the entire range 2 
T< B < 11 T. The REC and EXD ground-state energies 
exhibit oscillations as a function of B. These oscillations 
reflect the incompressibility of the many-body states as- 
sociated with magic angular momenta L m . The Z/ m 's are 
specified by the number of electrons on each polygonal 
ring; in general L m = Lq + J2q=i k q n q , with k q being 
nonnegative integers, i.e., L m = 6 + 4k2 for the (0,4) 
structure of N = 4, see Fig. 1. 

As a second accuracy test, we compare in TABLE I 
REC and EXD results for the interaction energies of the 
yrast band for N — 6 electrons in the lowest Landau 
level (an yrast state is the lowest energy state for a given 
L m ). The relative error is smaller than %0.3, and it 
decreases steadily for larger L values. The total energy 
of the REC is lower than that of the SEC, Fig. 1. A 
theorem discussed in Sec. Ill of Ref. ^£1, pertaining to 
the energies of projected wave functions, guarantees that 
this lowering of energy applies for all values of N and B. 

The crystalline polygonal-ring arrangement 
(n\, ri2, n r ) of classical point charges is portrayed 
directly in the electron density of the broken-symmetry 
SEC, since the latter consists of humps centered at the 
localization sites Zj's (one hump for each electron). 
In contrast, the REC has good angular momentum 
and thus its electron density is circularly uniform. To 
probe the crystalline character of the REC, we use the 



TABLE I: Comparison of yrast-band energies obtained from 
REC and EXD calculations for N = 6 electrons in the lowest 
Landau level, that is in the limit B — > 00. In this limit the 
external confinement can be neglected and only the interac- 
tion energy contributes to the yrast-band energies. Energies 
in units of e 2 /(kIb)- For the REC results, the (1,5) polygonal- 
ring arrangement was considered. For L < 140, see TABLE 
IV in Ref. [11(b)] and Table III in Ref. 0. 



L 


REC 


EXD 


Error (%) 


140 


1.6059 


1.6006 


0.33 


145 


1.5773 


1.5724 


0.31 


150 


1.5502 


1.5455 


0.30 


155 


1.5244 


1.5200 


0.29 


160 


1.4999 


1.4957 


0.28 


165 


1.4765 


1.4726 


0.27 


170 


1.4542 


1.4505 


0.26 


175 


1.4329 


1.4293 


0.25 


180 


1.4125 


1.4091 


0.24 


185 


1.3929 


1.3897 


0.23 


190 


1.3741 


1.3710 


0.23 


195 


1.3561 


1.3531 


0.22 


200 


1.3388 


1.3359 


0.21 
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allows calculation of the energies of RECs for arbitrary 
JV, given the corresponding equilibrium configuration of 
confined classical point charges. 

We focus on the description of the yrast band at a 
given B. Motivated by the aforementioned nonrigid 
character of the rotating electron crystallite, we con- 
sider the following kinetic-energy term appropriate for 
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configuration (with J2o=i ' 



JV): 



FIG. 2: CPDs for the ground state of N = 17 electrons at 
B = 10 T (L = 228). The observation point (solid dot) is 
placed (i) on the outer ring at ro = 1.858-Rq, left frame, and 
(ii) on the inner ring at ro = 0.969i?o, right frame. The rest 
of the parameters are the same as in Fig. 1. Lengths in units 
of Ro = (2e 2 / (nm* uiq)) 1 . CPDs (vertical axes) in arbitrary 
units. 



q=l 



h 2 L 2 /{2J q (a q )) - Hlu c L/2, (6) 



where L q is the partial angular momentum associated 



with the gth ring and the total L — ^2 a . 



1 Lq 



n q m*a q is the rotational inertia of each individual ring, 
i.e., the moment of inertia of n q classical point charges 
on the qth polygonal ring of radius a q . To obtain the 
total energy, Ef^ EC , we include also the term E^ pp ( JV) — 
Y^~ q =i Jq{ a q)Q 2 /2 due to the effective harmonic confine- 
P(r, ro) = ($| S(ri — r)S(rj — ro)|$) /($|$), (5) ment fi [see discussion of Eq. QJ], as wen as the interac- 



conditional probability distribution (CPD) defined as 



where $(ri, r2, rjv) denotes the many-body wave func- 
tion under consideration. P(r, tq) is proportional to the 
conditional probability of finding an electron at r, given 
that another electron is assumed at ro. This procedure 
subtracts the collective rotation in the laboratory frame 
of referenece, and, as a result, the CPDs reveal the struc- 
ture of the many body state in the intrinsic (rotating) 
reference frame. 

In Fig. 2, we display the CPD for the REC wave func- 
tion of N — 17 electrons. This case has a nontrivial three- 
ring structure (1,6,10) 0, which is sufficiently complex 
to allow generalizations for larger numbers of particles. 
The remarkable combined character (partly crystalline 
and partly fluid leading to a NCRI) of the REC is illus- 
trated in the CPDs of Fig. 2. Indeed, as the two CPDs 
(reflecting the choice of taking the observation point [ro 
in Eq. (J5J] on the outer, left frame, or the inner ring, 
right frame) reveal, the polygonal electron rings rotate 
independently of each other. Thus, e.g., to an observer 
located on the inner ring, the outer ring will appear as 
uniform, and vice versa. The wave functions obtained 
from exact diagonalization exhibit also the property of 
independently rotating rings (see e.g., the N = 9 case in 
Ref. which is a testimony to the ability of the REC 

wave function to capture the essential physics of a finite 
number of electrons in high B. 

In addition to the conditional probabilities, the 
solid/fluid character of the REC is revealed in its ex- 
cited rotational spectrum for a given B. From our micro- 
scopic calculations based on the wave function in Eq. J5J), 
we have derived (see below) an approximate (denoted as 
"app"), but analytic and parameter- free, expression [see 
Eq. (JSJ] which reflects directly the nonrigid (nonclassical) 
character of the REC for arbitrary size. This expression 



tion energy E^ pp , 



q=l 



EE** 

q=l s>q 



a.). (7) 



The first term is the Coulomb repulsion of n q point- 
like electrons on a given ring, with S q = X^=2( sm [(j — 
l)7r/n g ]) _1 , and the second term is the Coulomb re- 
pulsion between (uniform) rings with Vc(a q ,a s ) = 
n 9 n s2 Pi[3/4, 1/4; 1; 4a 2 q a 2 (a 2 + a 2 y 2 }e 2 (a 2 + a 2 s )-V 2 /K, 
where 2-Fi is the hypergeometric function. 

For large L (and/or B), the radii of the rings of the ro- 
tating crystallite can be found by neglecting the interac- 
tion term in the total approximate energy, thus minimiz- 



■E^ p (N). One finds a q 



^\fL q /i 



9) 



ing only E^(N) 
the ring radii depend on L q , reflecting the lack of ra- 
dial rigidity. Substitution into the above expressions for 
E^" p , E^ p and E^ yields for the total approximate 



energy the final expression: 



E 



REC 
app,L 



(JV) = K(n - co c /2)L + 



E 



c 



V,q 



-i T l/2 
q=l n Q 



EE W-' A \F)> (8) 

' ' \ n a V n s 



5=1 s>q 



where the constants Cy, q — 0.25n^ 2 S q e 2 / (k\) . In the 
case of the simplest (0, JV) and (1, JV — 1) ring configura- 
tions, Eq. JHJ reduces to the partial expressions reported 
earlier |14 l l20 | . 

In Fig. 3 (left frame), and for a sufficiently high mag- 
netic field (e.g., B = 100 T such that the Hilbert space 
of the system reduces to the lowest Landau level), we 
compare the approximate analytic energies E^ pp G L with 
the microscopic energies Ef^ EC calculated from Eq. © 
using the same parameters as in Fig. 1. The agreement 
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FIG. 3: Left: Approximate analytic expression [Eq. 
dashed line (online violet)] compared with microscopic REC 
calculations [Eq. solid line (online green)] for N — 17 
electrons at high magnetic field. Right: The corresponding 
classical (rigid rotor) energy £^ lg for N = 17 electrons (see 
text). The microscopic REC energies are referenced relative 
to the zero-point energy, 17?if2. Energies were calculated for 
magic angular momenta L = L\ + L% + £3 with L\ = 0, 
L 2 = 21 + 6fc 2 and k 2 = 30, and L 3 = 115 + 10fc 3 . The pa- 
rameters are the same as in Fig. 1 and B = 100 T. Note the 
much larger energy scale for the classical case (right frame), 
leading to a superfluidity index for the REC of ot a ~ 0.99 (see 
text). 



is very good (typically less than 0.5%). More important 
is the strong discrepancy between these results and the 
total energies E r ^ s of the classical (rigid rotor) crystal- 
lite, plotted in the right frame of Fig. 3; the latter are 
given by E 1 ^ = h 2 L 2 /(2J lig ) + 0.5 Zti m*w$\Zi\ 2 + 
'EriLi ^f>i e V ( K \Zi — Zj\), with the rigid moment of in- 
ertia being J lig = J2? =1 m*\Zi\ 2 

The disagreement between the REC and classical ener- 
gies is twofold: (i) The L dependence is different, and (ii) 
The REC energies are three orders of magnitude smaller 
than the classical ones. That is, the energy cost for the 
rotation of the REC is drastically smaller than for the 
classical rotation, thus exhibiting "supcrfluid" behavior. 
In analogy with Ref. we define a supefluidity index 
a s = (£^ s - Ef; EC )/E T i 6 . For the case displayed in Fig. 
3, we find that this index varies (for 1116 < L < 3716) 
from a s = 0.978 to a s = 0.998, indicating that the REC 
is highly superfluidic. 

Formation of a supersolid is often expected in conjunc- 
tion with the presence of (i) real defects and (ii) real 
vacancies 0, • Our supersolid wave function [Eq. [2] be- 
longs to a third possibility, namely to virtual defects and 
vacancies when the number of particles equals the num- 
ber of lattice sites (the possibility of a supersolid with 
equal number of particles and lattice sites is mentioned 
in Ref. [5}). Indeed, the azimuthal shift of the electrons 
by (71, 72j •••> 7r) [ see Eq. 1(2} ] may be viewed as gener- 
ating virtual defects and vacancies with respect to the 
original electron positions at (71 = 0,72 = 0, ...,7,- = 0). 

A recent publication |22J has explored the quantal na- 
ture of the 2D electron crystallites in the lowest Lan- 



dau level (B — > 00) using a modification of the second- 
quantized LLL form of the REC wave functions ■ ^ n 
particular, the modification consisted of a multiplication 
of the parameter free REC wave function by variation- 
ally adjustable Jastrow-factor vortices. Without consid- 
eration of the rotational properties of the modified wave 
function, the inherently quantal nature of the crystallite 
was attributed exclusively to the Jastrow factor. How- 
ever, as shown above, the original REC wave function 
[Eq. @] already exhibits the characteristic NCRI behav- 
ior of a supersolid. Consequently, the additional varia- 
tional freedom introduced by the Jastrow prefactor may 
well lead energetically to a slight numerical improvement, 
but it does not underlie the physics related to the super- 
solid behavior. Indeed, the essential step is the projection 
of the static electron crystallite onto a state with good 
total angular momentum [see Eqs. (JIJ and (J5}]. 

In summary, we have shown that rotating electron 
crystallites, formed in 2D quantum dots under high mag- 
netic fields, exhibit nonclassical rotational inertia which 
is the signature of a supersolid 0, Q . We have utilized 
an analytic many-body wave function [Eq. J5J)] which al- 
lowed us to carry out computations for a sufficiently large 
number of electrons (TV =17 electrons having a nontriv- 
ial three- ring polygonal structure), leading to the deriva- 
tion and validation of an analytic expression for the total 
energy of rotating electron crystallites of arbitrary N. 
These crystallites may be regarded as the precursors of 
an extended 2D fermionic supersolid in the lowest Lan- 
dau level. 
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